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The interplay between topological phases of matter and dissipative baths constitutes an emergent 
research topic with links to condensed matter, photonic crystals, cold atomic gases and quantum in¬ 
formation. While recent studies have suggested that dissipative baths can induce topological phases 
in intrinsically trivial quantum materials, the back action of topological invariants on dissipative 
baths has been overlooked. By exploring this back action for a centrosymmetric Dirac insulator 
coupled to phonons, we show that the linewidths of bulk optical phonons can reveal electronic band 
inversions. This result is the hrst known example where topological phases of an open quantum 
system may be detected by measuring the bulk properties of the surrounding environment. 

PACS numbers: 


Introduction.- The discovery of topological phases in 
three dimensional crystals has culminated in a new clas¬ 
sification scheme for solids that is based on quantum me¬ 
chanics and topology [ij. These phases are described 
by integers known as topological invariants, which mani¬ 
fest themselves through robust gapless states localized at 
sample boundaries. The characterization of topological 
invariants in insulators often idealizes electrons as being 
isolated from their environment. Yet, in real materials, 
electrons are coupled to various non-electronic baths and 
the usual idealization fails when the strength of the cou¬ 
pling exceeds the energy gap of the insulator. 

Recent work ii has suggested that baths can al¬ 
ter topological invariants and even induce topological 
phases. However, the inverse of this effect, concerning the 
back action of topological invariants on baths, remains 
completely unexplored. Does a change in the electronic 
topological invariant modify the surrounding bath? Is 
it possible to infer the topological invariants of an elec¬ 
tronic system by measuring its non-electronic environ¬ 
ment? The present work intends to answer these ques¬ 
tions affirmatively and thus establish an unanticipated 
interplay between band topology and dissipative baths. 

To that end, we adopt a minimal model which con¬ 
sists of massive 3D Dirac fermions coupled to a bath 
of phonons. In this model, we find that it is possible 
to learn whether the electronic band topology is triv¬ 
ial or nontrivial by analyzing the phonon linewidths in 
the thermodynamic limit (i.e. disregarding boundary ef¬ 
fects). Our results challenge a commonly held viewpoint, 
according to which the bulk properties of a doped topo¬ 
logical insulator and a doped trivial insulator should be 
qualitatively similar. 

Model. - The minimal Hamiltonian describing the low- 
energy bulk bands of a time- and inversion-symmetric 3D 
Dirac insulator near the Brillouin zone center is Q 

'H(k) = -f ak • crr“ -I- MkT^, (1) 


where cr* and r® are Pauli matrices in spin and orbital 
space (respectively), k = (kx,ky,kz) is the crystal mo¬ 
mentum, 7 models the particle-hole asymmetry of the 
band structure, a is the Dirac velocity, Mk = m + Pk'^ 
is the Dirac mass, 2|m| is the energy gap of the insu¬ 
lator at fc = 0, and P is an additional band parame¬ 
ter. Importantly, is the electronic parity operator and 
[t^= 0. For narrow-gap insulators described by 
Eq. dl]), the sign of m/3 determines the so-called strong 
topological invariant: if /3 > 0, then m > 0 (m < 0) re¬ 
sults in a trivial (topological) insulator. If m/3 < 0, the 
electronic bands at fc = 0 are said to be inverted. In addi¬ 
tion, Mk acts as a momentum-dependent effective mag¬ 
netic field that polarizes the orbital pseudospin t along 
the z direction. Because Mk changes sign as a function 
of k in the topological phase but not in the trivial phase, 
the fc-dependence of the expectation value of reflects 
the key difference between the bulk electronic structures 
of trivial and topological insulators (cf. Figs.[T^ and[Tl3). 

The eigenstates of Eq. © are |wkn), where n £ 
{1,..., 4} labels the two highest valence bands and the two 
lowest conduction bands near k = 0. The energy eigen¬ 
values, Ekn, are doubly degenerate owing to the com¬ 
bined time-reversal and inversion symmetries. For ana¬ 
lytical simplicity, we have chosen a continuum model with 
spherical symmetry; this captures the essential ideas and 
is smoothly deformable into more realistic lattice models 
that we use below for numerical calculations. 

Phonon self-energy.- Electron-phonon interactions 
shift phonon frequencies and contribute to phonon 
linewidth. These two effects can be calculated from the 
phonon self-energy Q, 


nA(q, WqA) 


_j_ (/krt /k-qra') /2\ 

V Ekn — Ek-qn' — WqA — f0+ 
knn ^ ^ 


Here, V is the sample volume, A labels different phonon 
modes, q is the phonon momentum, WqA is the bare 
phonon frequency and /kn is the fermion occupation 
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FIG. 1: (Color online) (a) and (b) Expectation value of the 
electronic parity operator, (r^) (represented by arrows), as 
a function of momentum for the electronic model of Eq. 0. 
(c) and (d) Fermi surface averages of |(t^)P and |{t“^)P (cf. 
Eq. 0), as a function of the Eermi energy, for m = —0.2 eV 
(c) and m = 0.2 eV (d). The rest of the band parameters are 
the same as those in Ref. [^. 

number for the state |ukn) with a Fermi energy ep. Also, 
5L'(k,q) = (wknlff^(q)lwk-qn'), (3) 

where g^{q) = ff'''(—q)^ is the electron-phonon vertex 
operator in the low-energy electronic subspace [1, . 

In a centrosymmetric crystal, lattice vibrations are ei¬ 
ther even or odd under spatial inversion. Each of these 
modes couples to electrons and can, as we shall see, in¬ 
herit signatures of the underlying band topology. For the 
model of Eq. © and for g ~ 0 optical phonons, inversion 
and time-reversal symmetries dictate [g 

^even(q) _ gg(q) -p 5^(q)r^ 

g°'i'i(q):^g,(q)r"+g'(q).<Tr^ (4) 

where “even” (“odd”) denotes the coupling of elec¬ 
trons to parity-even (parity-odd) phonon modes, with 
_ q _ q_ Inversion symmetry 

guarantees that and g°‘^^ will not be mixed in a sin¬ 
gle phonon mode. Also, q = q/g, and the coefficients gi 
{i = 0, X, z) and g' (f = x, y, z) can be obtained from the 
atomic displacements in the particular phonon mode @ . 
Physically, go and g^ lead to phonon-induced modula¬ 
tions of the chemical potential and the Dirac mass, re¬ 


spectively. Next, we identify ways in which g^ can trans¬ 
fer the information about electronic band topology to the 
phonon sector. 

Intraband phonon damping.- The main electronic 
mechanism contributing to phonon linewidths is the scat¬ 
tering of phonons off electron-hole pairs. The rate of 
this process is 7^(q) = —ImnA(q,WqA). In this work, 
we focus on long-wavelength optical phonons and on low 
temperatures. 

We begin by considering the commonly realized case in 
which the phonon frequency is smaller than the bandgap 
of the insulator. In this case, the “insulator” must be 
doped in order for carriers to absorb phonons and in¬ 
duce a linewidth The subscript “ia” is shorthand for 
“intraband” and makes it explicit that phonons decay 
into particle-hole pairs in the vicinity of the Eermi sur¬ 
face. Assuming that the distance from the Eermi level 
to the bulk band et^e is large compared to the phonon 
frequency, we have [g 

7il(9 ^ 0) ~ TruJo\D{€F)\g^^{kFAWS{vF • q - wqa), 

( 5 ) 

where D{eF) is the electronic density of states per band 
at the Fermi level, kj’ is the Fermi wave vector, -vp is the 
Fermi velocity, 6{x) is the Dirac delta, and q)P de¬ 
notes the sum of over the two degenerate bands 

at momentum k and energy Ak (hence the label “intra¬ 
band”). In addition, O = ~ ^F)/{VD{eF))- 

Equation © contains information about the electronic 
band topology. The simplest way to see this is to imag¬ 
ine a parity-even phonon mode and a parity-odd phonon 
mode that couple to electrons purely through g^ = g^r^ 
and g^ = gxT^, respectively. More general couplings with 
go ^0 and gi ^0 will be discussed below. Erom Eqs. o, 
© and the linewidths of these two phonon modes 
are 0 

7 , 4(9 ^ 0 ) ~ |gj(q)pD(eF)|('rJ)iap Y ®(1 “ (®) 

where j € {x,z}, 0(x) is the Heaviside function, g = 
i^oj/iqvF) and 

TPhTP = 1 - (7) 

Note that |(r-?)iap G [0,1] (cf. Fig.©. In particular, 
when Mkp = 0, |(T^)iap = 0 and |(T"')iaP = 1. Com¬ 
bining Eqs. m and © with Eig. [II it follows that 
reflects the orbital texture, and therefore the topology, 
of the bulk bands. In order to clarify this point, we elim¬ 
inate the non-topological features coming from D{eF) by 
considering the ratio - (9^/5'Dl('r"^)ia| Vl('r^)iap- 

Eor fixed bandgap, Eq. © predicts a strong maximum 
for yia/yil as a function of ep in the topological phase 
(but not in the trivial phase) because Mkp crosses zero as 
a function of ep in the topological phase (but not in the 
trivial phase). This difference in behavior between the 
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trivial and topological phases is significant for a sizeable 
|m|, but becomes gradually weaker as the energy gap 
decreases, ultimately disappearing when m —>■ 0. In other 
words, 7 ia/ 7 fa contains no signatures of band topology 
near the topological quantum critical point. Figure [2^ 
confirms our analytical statements in a lattice model, for 
which Eq. is solved numerically. 

Alternatively, in a sample with fixed carrier density, 
7 ia/ 7 fa shows a pronounced maximum as a function of m 
in the topological phase only. The maximum takes place 
at m* ~ —Pkp, where Mkp undergoes a sign change. 
Motivated by recent claims of pressure-induced band in¬ 
versions in Sb 2 Se 3 and Pbi_a;Sna;Se [ll|) [s3|) in Fig. [Hd 
we plot 7 ia/ 7 fa as a function of pressure, using a lattice 
model. This corroborates the emergence of a “topology- 
induced” maximum in 7 ? /yf,. 

In the preceding discussion of 7 ;“, we have assumed a 
parity-odd phonon mode that couples to electrons purely 
through {g[ = 0 in Eq. (jH)). In general, such a phonon 
can also couple to electrons through the term g' • crr^. 
However, we have verified that this coupling produces 
qualitatively similar features as r^. 

Similarly, when discussing we have imagined a 
parity-even phonon mode that couples to electrons purely 
through {go = 0 in Eq. ([5])). Nonetheless, symm etry 
allows a mixture of and the identity matrix 1 [l3j |. 
The latter produces intraband matrix elements that are 
insensitive to the orbital texture of the insulator, since 
(Mkn|l|ukn) = 1- Consequently, the effect of 50 7 ^ 0 is to 
dilute away the topological features of 7 ?^. Although this 
constitutes a problem towards the realization of Fig. [2] 
in real materials, we find that the maximum in 7 ia/ 7 ia 
remains pinned to the topological side if \gz\ > |5o|- 

Interband phonon damping. - Thus far, we have consid¬ 
ered the linewidths of phonons with woa < 2|m|. Herein, 
we investigate the case wqa > 2 | 77 i|, relevant to Dirac 
insulators with particularly small bandgaps and/or high- 
frequency phonon modes. In this case, a phonon is ab¬ 
sorbed by an electron in the bulk valence band, which 
gets promoted to the bulk conduction band. The asso¬ 
ciated phonon linewidth is 71 ^, where the subscript “ie” 
is shorthand for “interband”. Assuming for the moment 
that cf is inside the bulk gap, Eq. (1551) yields [1] 

7i^((7 ~ 0 ) ~ 7rDjoint(woA)|ffie(k,q)P, ( 8 ) 

where Djoint(t^) = Sk'^(^kc - E]^v - is the joint 

density of states, Akc and are the bulk conduction 
(c) and valence (v) band energies. In addition, IfteP = 

Enec,n'ev l5wP and \g{'J^ = Ek - Ak„ - 

Woa)/ (VDjoint)- 

From Eqs. o and we obtain |( 7 ;g(k, 0 )p = 

l 5 fa(k, 0 )P and |gfg(k,0)|2 = |gf^(k,0)p. Therefore, 7 ;^ is 
as sensitive as 7 ^ to the band topology of the Dirac insu¬ 
lator (with X and 2 ; interchanged). More so, an important 
advantage of 7 jg over is that we may effectively take 
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FIG. 2: (Color online) (a) y^i/yfa as a function of the Fermi 
energy and the bulk carrier density, for m = ±0.1 eV. A 
prominent maximum emerges in the topological phase only, 
due to the momentum-space orbital texture of the electronic 
eigenstates, (b) 7 ^ 1/74 as a function of pressure P. We use 
m = a{P — Pc), where Pc is the critical pressure for a band 
inversion and a is a coefficient that can be obtained e.g. from 
experiment fs^ . The bulk carrier density is n ~ no(l±P/B), 
where no is the density at P = 0 and B is the bulk modu¬ 
lus. The maximum of yjl/yfa, appears at P = P* . Inset: The 
dependence of P* on no. As no decreases, P* approaches 
Pc, making it more difficult to identify trivial and topolog¬ 
ical phases solely from phonon measurements. Throughout 
this figure, we have used a tetragonal lattice regularization 
of Eq. ([ 1 ]). Because q;,/ 3,7 are not tabulated for Sb 2 Se 3 , we 
have replaced them with those of SbjTes 0]. For the bulk 
modulus, we have used B = 30GPa [®| . 


-even _ -z of the value of go in Eq. ®, because 

(Mkn|l|Mkn') = 0 for iutei'band transitions. Accordingly, 
the topological signatures in yj^ are more robust than 
those in y^- 


In a sample with fixed carrier density, 7ie/7ie contains 
a minimum as a function of m at m* ~ —a;g/3/(4a^), 
i.e. only in the topological phase 1^. This result has 
the same origin as the maximum of 7ia/7fa discussed 
above, and it holds for doped samples as well so long as 
akp/ujo -C 1. Figure [5] confirms this for a lattice model. 


Discussion. - In sum, there are three reasons why the 
linewidths of bulk, long-wavelength optical phonons can 
inherit distinct signatures of the electronic band topol¬ 
ogy. First, phonon linewidths are proportional to the 
square of energy-resolved electron-phonon matrix ele¬ 
ments. Second, in a centrosymmetric crystal, the cou¬ 
pling of a optical g ~ 0 phonon to electrons either com¬ 
mutes or anticommutes with the electronic parity oper¬ 
ator. Third, the momentum-space texture of electronic 
parity eigenvalues reflects band inversions. 


To be precise, the predicted features in the phonon 
linewidths probe band inversions near the Fermi level, 
rather than the strong topological invariant per se. How¬ 
ever, in materials whose low-energy bands are described 
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FIG. 3: (Color online) 7ie/7ie S’® a function of the Dirac 
mass m, where the rest of the band parameters correspond to 
Sb 2 Te 3 (a) or Bi 2 Se 3 (b). The minimum of 7ie/7fe occuring 
in the topological side is a direct manifestation of the orbital 
texture in Fig. [T] Throughout this figure, we have used a 
tetragonal lattice regularization of Eq. m with the band pa¬ 
rameters taken from Ref. @]. The Fermi level is assumed to 
be inside the bulk gap. 


by Eq. O, the energy gap minimum occurs at a sirrgle 
time-reversal-invariant momentum (TRIM) and hence a 
band inversion taking place therein is equivalent to a 
change in the strong topological invariant. By exten¬ 
sion, phonon linewidths can probe the strong topological 
invariant in centrosymmetric crystals where the direct 
bandgap minimum is known to occur at an odd num¬ 
ber of symmetry-equivalent TRIM. In contrast, phonon 
linewidths do not faithfully reflect the strong topological 
iirvariant iir materials where the direct bandgap irriirimum 
can occur at an even number of symmetry-equivalent 
TRIM, nor in materials lacking an iirversion cerrter. 

Phonon frequencies, which we have barely mentioned 
thus far, are much less sensitive than phonon linewidths 
to the electronic band topology. This is because the real 
part of Eq. (1331) contains a sum over electron-phonon 
matrix elements at multiple energies, with weights that 
depend on non-topological details of the energy bands. 
This notwithstanding, a recent experiment in Sb 2 Se 3 
has attributed a kink in the pressure-dependence of the 
phonon frequency to a band inversion. Our calcula¬ 
tions @ do not support such interpretation. 

The main tools to measure g ~ 0 phonon linewidths are 
Raman spectroscopy (for parity-even phonons), infrared 
spectroscopy (for parity-odd phonons) and inelastic neu¬ 
tron scattering [l5|. In a clean material with u!o\t ^ I 
(where t is the disorder scattering time), vanishes 
unless q > ujo\/vf [I^. Since ujqx/vf typically exceeds 
the photon wave vector used in optical spectroscopies, 
7 ^ should be measured with neutrons. In contrast, 7 ;g 
remains nonzero at g = 0 and is thus amenable to optics. 
For Bi 2 Se 3 , we estimate ^ la ^ 1 cm“^, which nears the 
experimental resolution [l7j |. 

Aside from electron-phonon interactions, anharmonic 


lattice effects contribute to the phonon linewidth. To 
leading order, phonon-phonon interactions contain no in¬ 
formation about the electronic band topology and are 
independent from the carrier density. Therefore, the 
anharmonic part can be subtracted by measuring the 
linewidths with respect to a baseline carrier density. 

In view of our results, it is natural to ask whether 
any other physical observable involving Fermi’s golden 
rule, such as conductivity, might be sensitive to elec¬ 
tronic band topology on the same footing as the phonon 
linewidths. The answer is generally negative. For exam¬ 
ple, the optical conductivity cannot clearly differentiate 
between trivial and nontrivial orbital textures in the bulk 
because the velocity operator mixes 1, and [l8| . 

To conclude, we have proven that it is in principle pos¬ 
sible to infer the strong topological invariant of multi¬ 
ple Dirac insulators from the linewidths of bulk, long- 
wavelength optical phonons. It may be of interest to 
investigate formal links between the phonon linewidth 
and the SU(2) Berry phase identified in Ref. 0. Like¬ 
wise, it will be desirable to complement our theory with 
ab-initio electronic structure calculations, and to search 
for similar insights in other contexts (cold atoms, pho¬ 
tonic crystals, quantum memories) where the interplay 
between topology and dissipation may be crucial. 
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SUPPLEMENTARY MATERIAL 


Electron-phonon conpling 

The objective of this section is to derive Eq. (4) of the main text. For the first part of the derivation, we follow 
Ref. [sH . We start with the periodic lattice potential given by 


C/(r)=^t/^(r-R(l,/3)), 

1./3 


( 9 ) 


where R(/,/3) = \ + Tp is the equilibrium position of an atom (3 located in the th unit cell, 1 being the position of 
the unit cell and Tp being the position of the atom with respect to a reference point in that unit cell. The departures 
from the atoms from their equilibrium positions produce a change in the lattice potential. 


SU{r) = Y, Up{v - R(l, p) - Q(l, P)) -YUp{v- R(l, P)) 

1./3 1./3 

dUp{y-^{\,P)) 


( 10 ) 


where Q(l,/3) is the displacement operator for atom P in unit cell 1 away from its equilibrium position, given by 

= E . e--^-9+-^)pA(q,/3)(a-qA + 


q,A 


2MpNuJq\ 


( 11 ) 


Here, q is the phonon momentum, aqx is an operator that annihilates a phonon mode A with momentum q, Mp is 
the mass of atom /3, PA(q,/?) is the polarization vector corresponding to atom P in mode A and momentum q, WqA is 
the phonon frequency and N is the number of unit cells in the crystal. 

Then, the local electron-phonon coupling can be written as 


= J drp{r)SU{r), 

where p(r) = '0l(r)i/)(r) is the electron density operator, 

V’(i’) = E e“*‘''‘'^kn(r)ckn, 


( 12 ) 


(13) 


V is the volume of the sample, Ukn(r) is the periodic part of the Bloch wave function corresponding to band n and 
momentum k (within the first Brillouin zone), and Ckn is an operator that annihilates an electron in such a state. It 
follows that 


1 




J(k-k').r * 


Kn(r)uk'n' (r)Q(l,P) 


dUp(r-R(l,P)) t 


/31 kk' nn' 


dR(i,p) 


^kn^k'- 


=^EEE / rfre*('^-‘^') (--+')4Jr + l)uk'n'(r + l)Q(l,/?)- 


9U;3(r + l-R(l,/3)) ^ 


N 


/91 kk' nn' 


dR(i,p) 


C-kn^k'n' 


=vEEEE^- 


k+q1 


kk' nn' qA 


■EEEE 


Vcell ^ ^ ^ V 2MpNu}qX 

k nn' qA y 


J dreg'll '■+'^''^wL(i')uk'n'(r)pA(q,/3) • E^E^^^EMct^Ck'„'(a-qA+ aqA) 

//3(r-R(0,/3)) ^ ^ ^ I ^ 

aR(0,/3) '^knCk+qnKa-qA + flqA) 


2MpNujqX 

J dre*'i (”'’+'^'5^Uk„(r)uk+qn' (r) PA(q, P) 


9nn' (kj q)Ck„Ck+qn'(®-qA + Qq^)' 

^ knn' qA 


(14) 
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where Vceii = V/N is the unit cell volume and 


9nn 




2M/3WqA 


J ( '■+'^'5^<„(r)uk+q„/(r)pA(q,/3) 


dU0ir-R{O,(3)) 

9R(0,/3) 


Vv< 


cell 


E 


h 


‘^MpLUqX 




aR(0,/3) 


(15) 


is the matrix element that appears in Eqs. (2) and (3) of the main text. In Eq. (fni) . we have used u{r +1) = M(r) and 
exp(ik • 1 ) = iV(5k,o- It is implicit that, whenever k + q exits the first Brillouin zone, a reciprocal lattice vector will 
be added to bring it back in. The units of are energy x Vvolume. In spite of the integration over the entire crystal 
in Eq. 113, is an intensive quantity (i.e. independent of the volume of the crystal in the thermodynamic limit). 
This is a consequence of the fact that 9?7/3(r — R(0,/3))/9R(0, /?) is localized in space, while the Bloch functions are 
extended. 

Thus far everything has been general |Sll |. From now on, we will limit ourselves to the low-energy electronic bands 
spanned by the two highest valence bands and the two lowest conduction bands in the vicinity of the E point. Within 
this four-band subspace, it is useful to transform 5 ^„,(k, q) into the basis spanned by lucrr), which are the low-energy 
electronic eigenstates at T (i.e. at fc = 0). Specifically, 

l^kn) — ^ ) I^crr) (^crr I'akn) ; (^^) 

<TT 

where a and r are the two-level degrees of freedom describing spin and orbital, respectively, and (uo-rlukra) can be 
obtained by diagonalizing a k • p Hamiltonian (e.g. Eq. (I) from the main text). Then, 

ffnn'(^!*l) ~ 'y ^'^~~! ('^kra|Mo-T)('a<T'T'|^k+qra')go-T:cr'T'( q)) (1'^) 

err cr'r' 


where 








cell 


2M/3UJq\ 


dre *'i-'’u<^r(i’)*w<T'T'(r) PA(q,/3) 


dU0{r-R{O,/3)) 
5R(0,/3) ■ 


(18) 


It is now possible to rewrite in the basis spanned by cr and r: 

<^=4eeee5^^- '(q) CkCTr''k-(-qcr'r'(a-qA + OqA)' (1^) 

^ k qA err a'r' 

In general, the evaluation of g^^.^,^,{q) requires a detailed knowledge of the crystal structure of the material and of 
the electronic wave functions at F. Here, we are interested in making some generic statements based on symmetries. 
Let g^(q) be a 4 X 4 matrix such that 

{UaT\g^{q)\Ua’T') = 3A;a'r'(q)- (20) 

When the matrix g^(q) acts on |ucrr)j the latter can be thought of as a 4-component spinor. The most general form 
for g^(q) is 


ff^(q) = ^ 

u 


( 21 ) 


where i,j £ {0,x,y, z}, the 0-th Pauli matrix represents the identity, and 5 )^ (q) are complex numbers. There are 
various constraints on the form of the coefficients g^y First, since TL^gp is Hermitian, the matrix g must obey 

5"(q)=ff"(-q)^ (22) 


Second, time-reversal symmetry dictates 


0-V(q)0 = p^(-q), 


( 23 ) 
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where 0 = ia^K is the time-reversal operator and K is the complex-conjugate operator. By combining Eqs. (1221) and 
(l23l) . it follows immediately that 

only matrices that commute with 0 are allowed in the expansion of g. (24) 

This rules out many terms in Eq. (ED). Yet, as we show next, if the crystal has inversion symmetry, there are additional 
constraints in the limit of long-wavelength phonons {q — 0), which simplify the form of g considerably. When g ~ 0, 
Eq. (1181) may be simplified through 

~ 1 (25) 

and 


~ 1 — iq • r 


(26) 


In Eq. (1^51) . the Taylor expansion has been truncated at the leading term because the size of a unit cell is negligible 
compared to the phonon wavelength when g ~ 0. The expansion in Eq. (E51) is at first glance less straightforward, as r 
can be arbitrarily large. However, let us recognize the fact that 9?7(r —R(0, ,8))/i9R(0, /?) is a function that is localized 
about the 1 = 0 unit cell. Then, the main contribution to the integral in Eq. (TT^ comes from r in the “vicinity” of 
the 1 = 0 unit cell, where “vicinity” is quantihed by the range of dU{Y — R(0,/3))/i9R(0,/3). For g ~ 0 phonons, this 
range is always short compared to the phonon wavelength, which then justifies the expansion in Eq. (1261) . Moreover, 
for <7 ~ 0 optical phonons, it is sufficient to keep only the leading term in Eq. (1261) . A similar type of argument has 
been invoked in Ref. S2|. 

The main advantage of considering the long-wavelength limit in a centrosymmetric crystal is the emergence of a 
simple parity selection rule. Let us begin by writing Eq. dUD as 


9 ar 


;a'r'(q) -X! 



(r)pA(q,/3) 


dUp{v-^{0,p)) 


'(-i')Pa(4/3) 


aR(0,/3) 

dUp{-v -^{0,13)) 

5R(0,^) 


(27) 


Now, we exploit the fact that \u^t) is an eigenstate of the parity operator: = {—^YWaT)- Choosing r = 0 at 

the center of inversion, we have 


9 ( 7 T\a 


,,(q) =(-!)-+”'^ 
P 

=(-ir"'E 


-(-r) = (-l)”'u„^(r). Then, 

1 I h 


A/Vcell y 

1 I h 


dr uY{r)u„>^>{r)pxi(i,/3) 


dUp{-r-R{0,p)) 


■\/Vceii V 2M_/3a;qA 


J druY{r)ua>r’{r)pxiq,-/3) ■ 


5R(0,/3) 

dU.p{-r-R{0,-l3)) 
aR(0, -/3) 


(28) 


where, in the second line, we have simply changed the dummy index /3 to — /3. Now, we recognize that (i) /3 and 
—/3 are the same atom (related by the inversion operator), (ii) R(0,/3) = —R(0, —/3) (because we take the origin 
at the center of inversion), (iii) Upi^—r R(0,/3)) = Up{v — R(0,/?)) (because Up{Y — R(l,/3)) = Up{\v — R(l,/3)|), 
(iv) PA(q, —/3) = (—l)^~''^PA(q,/3), where we set the convention that A = 0mod2 for parity-even phonon modes and 
A = lmod2 for parity-odd phonon modes. This last relation emerges from the fact that, in a g = 0 parity-even 
(parity-odd) mode, two atoms related by inversion have opposite (equal) displacement vectors ( see Ref. [S3j | ). Then, 


9ar:a'T' 


(q) =(-i) 


r+r' 


= (-l 




T 

p 

T 


1 


\/Vcell y ^Af—^CJqA 

1 I h 


J rfr<r(l’)M<T'r'(l’)PA(q,-/?) 


\/Vcell y 2M,gWqA 
1 


J druY{r)ua>r'Y)Px{q,-P) 


=(-i)”+”'(-i)^+i(-i)E 


h 


y/p-f, y 2M^WqA 


druY{r)ucr'T'ir)px{q,l3) 


dU.pi-r-R{0,-P)) 
5R(0, -(3) 

dUp{r-R{0,(3)) 

(-)aR(0,/3) 

dUp{r-R{0,f3)) 


dR{0,13) 


=(-l) 


r+T +A A 


9gT-,(7't' (^)* 


(29) 
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Thus, the parity selection rule reads t + t' +A = 0mod2. In words, a parity-even g ~ 0 phonon mode can only couple 
electronic states of the same parity, while a parity-odd g ~ 0 phonon can only couple electronic states of opposite 
parity. These statements are no longer exact when we keep the iq- r term in Eq. (1261) : however, the correction due to 
such term is necessarily small when g ~ 0. 

The parity selection rule from Eq. dMl), together with Eq. (1211) . implies that 

= Soo(q) + 9oz{^)t^ 

r‘'‘'(q)=50x(q)r"+ ^ g.,(q)aV^. (30) 


Namely, the electronic operators appearing with parity-even phonon modes commute with t'^, while those appearing 
with parity-odd phonon modes anticommute with . This is the most important result from the Supplementary 
Material: it constitutes Eq. (4) of the main text (N.B. in the main text we have lightened the notation slightly. 
For example, we have used go rather than goo, Qz rather than goz, etc.). The coefficients appearing in Eqs. (|3nil 
can be computed from first principles with the aid of Eq. dUl). For instance, goo = (1/4) Ear 5oz = 

(1/4) Ear(-l)"5a^;7r", 30r = (1/4) Ear 5ar: etC. 

Up until now we have focused on a spatially local electron-phonon coupling. Let us briefly comment on the nonlocal 
electron-phonon coupling. To leading order, we may model the non-local electron-phonon coupling by 

^non-loc dv5U{v)i,\v)^{v + d), (31) 

<5 

where 8 is the vector connecting nearest neighbor atoms. Following the same steps outlined for we obtain 

= ^5^^^gE,(k,q)e-‘(‘^+'i)'^cLck+q„'(a-qA + a^,). (32) 

^ kqA nn' 5 

Eqs. (fT4l) and (1321) are formally similar, except for the factor e“4k+q)'5 latter. An immediate consequence of this 

extra factor is that g^ becomes a function of both q and k rather than just q, and that there are more allowed terms 
than those written in Eq. ((30)l . Nevertheless, we are interested in weakly doped Dirac insulators, which means that 
the momentum of the electrons contributing to the phonon linewidths is small compared to the size of the Brillouin 
zone. This, together with the fact that 8 is comparable to the inverse of the size of the Brillouin zone, implies that 
k-^ ^ 1. Thus, the nonlocal electron-phonon interaction approximately reduces to the local one and Eq. (1301) remains 
a good approximation for long-wavelength phonons. 


Real part of the phonon self-energy 

In the weak coupling regime, the electronic-induced changes in phonon frequencies are proportional to the real part 
of the phonon self-energy. 


RenA(q,WqA) 


_J_ q)P(/k" /k-qri') 

V Ekn T/k—qn^ ^qA 

k,n,n 


(33) 


Equation (1331) involves a sum over electron-phonon matrix elements (which contain information about the orbital 
character, and hence topology, of the Bloch bands) weighted with a function that depends on the details of the energy 
dispersion (which is only marginally related to the topological invariant). It follows that the topological features of 
RoIIa are rather diluted. This is in stark contrast with ImllA, studied in the main text, where the presence of the 
Dirac delta gives way to energy-resolved electron-phonon matrix elements. 

In support of the preceding observations. Fig. illustrates RenA(0,a;oA) as a function of bare Dirac mass m for 
the “a; and 2 modes” discussed in the main text (which couple to electrons via and r^, respectively). In addition, 
for reference purposes, we plot RenA(0,a;qA) for a constant matrix element (gE'(k, q) = 1). The similar qualitative 
behavior for the “z mode” and for the case of constant matrix elements corroborates our view that RellA is weakly 
sensitive to the orbital texture of the band structure. We also notice that the behavior of the “x mode” is somewhat 
different. This difference does originate from the electron-phonon matrix elements and does bear somewhat on the 
orbital character of the eigenstates; unfortunately, there are no prominent model-independent features that would 




FIG. 4: (Color online) Real part of the electron self-energy for q = 0 optical phonons. The red (dot-dashed) and the blue 
(dashed) curves correspond to the “x” and “z” modes, respectively. The black (solid) line corresponds to the case where 
p„„/(k,q) = 1. The band parameters (Sb 2 Te 3 ) are borrowed from the Ref. [4] of the main text. Also, a;oA = lOmeV and 
ep = 50 meV. 


allow us to learn about the band topology of a Dirac insulator by comparing the frequencies of the x modes and the 
z modes. 

We conclude by stressing that RcIIa is smooth across a band inversion in doped samples. This statement is in 
apparent disagreement with Ref. |S4j . where it has been claimed that RcIIa displays sharp features or anomalies at 
(or close to) a band inversion. Moreover, we find that the non-monotonicities in RellA (e.g. its changes in slope as a 
function of m) are model-dependent and not particularly reflective of the underlying band topology. 
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Imaginary part of the phonon self-energy 

In this section, we derive Eqs. (5), ( 6 ) and ( 8 ) of the main text. The starting point is the imaginary part of the 
phonon self-energy, given by 

ImnA(q,a;qA) ^ X! l5w(k, q)P(/kn - /k-q«')'^(^kn - Ek-q„' - WqA). (34) 

knn' 

We limit ourselves to (i) long wavelength optical phonons, (ii) a four-band electronic structure with pairwise degenerate 
conduction and valence bands, and (iii) low temperatures. 

Let us first suppose that the phonon frequency is smaller than the bulk energy gap. In this case, 5{E\^n — Ek-qn' ~ 
WqA) vanishes unless n and n' belong to energy-degenerate bands. Consequently, 

X! l5nn'(k,q)P 

knn^ 

^kn—^'kn^ 

-l 5 w(k,q)p( 5 (vkn • q-W 0 A)( 5 (Ek„ - cf), (35) 

knn^ 

Ekn. = Ei^„l 

where Vkn = 5ifkra/9k. The expansion of the Fermi function in the first line of Eq. (1551) is valid as long as the distance 
between the Fermi level and bulk band edge is large compared to the phonon frequency Wqa- Note that g > 0 is needed 
in order to have a non-vanishing rhs of Eq. (I35p . 

The two bands that intersect the Fermi level are degenerate (due to inversion and time-reversal symmetry). Let us 
define their energy dispersion to be E\^. It follows that Vkn = dE^/dh = Vk for both values of n. Then, we may write 

ImnA(q ~ 0 ,a;qA) ^ |g 4 (kF, q)p5(vF • q-a;oA)^(Ek - ep) 

k 

= - 7 rwoAL>(eF)|gii(kF,q)p( 5 (vF • q- wqa), (36) 

where vf = Vk^. is the Fermi velocity, D{€f) = p X(k ^(-^k — cf) is the density of states per band at the Fermi level, 

lffia(kF,q)P = Y l5w(k,q)P (with Fk„ = Fk„- = cf), (37) 

nn' 


TT 

ImnA(q ~ 0 ,a;qA) ^ — 


/kn /kn T ^OA 


dhn 

dEi^ri 


(5(vkn -q-WoA) 


and 


0 = ^0k5(Fk-eF)/^5(Fk-eF) (38) 

k k 

is the Fermi-surface average of O. Hence, Eq. (5) in the main text is obtained via the relation 7 ^ = 
-ImnA(q ~ 0,a;qA). 

From Eq. (IM)) . it is straightforward to derive Eq. (6) of the main text. For this derivation, we adopt the model 
Hamiltonian of Eq. (I) in the main text, which is spherically symmetric. Accordingly, we may write 

ImnA(q~ 0,a;qA) X! Is'ialk, q)P<5(vF -q - —)(5(Fk - cf) 

VqvF QVp 

k 

poo /*1 p27r 

= - — pTAa / d{cos9) d(p \g^^{k,ci)\^ 6 {cos 9 - r]) 6{Ei, - ep), (39) 

qvp Jo {2 ttY j_i Jo 

where Vp = Vp/vp, q = q/g and rj = wq\/{ qvp). At this point, we are interested in phonon modes that couple to 
electrons through r® or (cf. the paragraph preceding Eq. (6) in the main text). For these, we can write 

lffia(k,q^ 0)p = |c,j (q)|2 ^ |(Mkn|T^ |Mkn')|^ = |gj(q)n(T^)iaP, (40) 

nn 

^kn = ^kn' 
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where the sum over n and n' is constrained to the two degenerate bands that intersect the Fermi surface, and j 
equals either x ot z (N.B. there is no sum over j). In the spherically symmetric model, |(T'^)iaP depends only on the 
magnitude of k. Consequently, 

Imnj-(q cTi; 0,a;qj) = 27r / - ep) f d{cos9) S {cos 9 - r]) 

qvF Jo (27r)-^ 

= _ 7ra;oj|gj(q)| f d{cos 9) 6 {cos 9 - ij) 

QUf Jo {ztt) 7-1 

= _ 7ra;oj|gi(q)| 27r|(Tt)iaP-|^ f d{cos 9) S {cos 9 - rj) 

_ 97ri-2 

=-7rr;|g,(q)n(r^)iaP—^0(1 - g) 

[Ztt 

= -7rg|gi(q)n(T^')iaP ^^^^^ 0(l - g), (41) 

where 0(a:) is the step function and 


-P(e-F) = V ~ = ir, yi —■ (42) 

V (27r)'^ VF 

k 

By taking the negative of Eq. (41), we obtain Eq. (6) of the main text. 

Finally, we derive Eq. (8) of the main text. In order to do so, we consider the case where the phonon frequency 
exceeds the bulk energy gap. Then, 

ImnA(q ~ 0,WqA) ^ |g^„,(k,q)p(/k„ -/k„/)^(£'k„ - Ekn'- woa), (43) 

k,nGc,n'G'i’ 


where c and v correspond to conduction and valence bands, respectively. Defining the respective energies as Eke and 
Eki,, and assuming for simplicity that the Eermi energy lies inside the Dirac gap, we have 


ImnA(q ~ 0,WqA) l5ce(k, q)p(5(Ekc - Ek^ - wqa) 

k 

= -7rEjoi„t(a;oA)|gie(k,q)P, 


where IsL'P. 

.^joint(a^OA) — ^ ) ^(Ekc Eke, ^Ioa) 

k 

is the joint density of states and 

l5i4(k,q)P = Y l5ie(k, q)P^(-E'kc - Ek„ - a;oA)/X!'^(^kc - Ek„ - wqa)- 

k k 


(44) 


(45) 


(46) 


Thus, we recover Eq. (8) in the main text. 
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